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Abstract 

We discuss the resummation of Sudakov effects in a cross section from the 
viewpoint of its underlying factorization near the edge of phase space. We 
perform the resummation of Sudakov threshold logarithms in Higgs (Standard 
Model and Minimal Supersymmetric SM) production at the LHC, using an 
evolution equation in the Higgs mass that is derived from this factorization. 
We extend the class of universal large terms that is resummed to include ad- 
ditional universal contributions, which, when included, help to reproduce the 
exact result to within a few percent for the full allowed range of Higgs bo- 
son masses in the SM and MSSM. Using the analytic resummed formula as a 
generating functional for approximate perturbation theory, we show results for 
next-to-next-to-leading order corrections in Higgs production, and find they are 
potentially sizable. 



1. Sudakov Factorization and Resummation 

The intimate connection between resummation of large logarithms in amplitudes 
in quantum field theory, and the factorization of such amplitudes is well-known from 
multiplicative renormalization. The unrenormalized Green function in terms of fields 
4>i factorizes as 

G nn (pi, A, g Q ) = Yl \[Zi(-, g(n)) G rcn (pi, /i, g(/i)) (1) 

i ft 

where A is a UV cutoff, and go/g{f^) is the unrenormalized/renormalized coupling. 
G un does not depend on /x, nor does G rcn on A, so that 

^ = -E7^)) (2) 

where ji = (/xd/ dfx) ha. y/Zi can only depend on the renormalized coupling, by sepa- 
ration of variables. The solution to this evolution equation (here the renormalization 
group equation), is the resummed Green function. In this example single logarithms 
are resummed. We wish to exhibit a similar paradigm for the double logarithmic, 
or Sudakov, case. The connection between factorization and Sudakov resummation 
was already pointed out by Mueller, Collins and Sen [TJ. We only provide a brief 
description here of a streamlined approach JJ] to Sudakov resummation for various 
reactions, that stresses common features. 

We consider cross sections near the edge of phase space, where there is not much 
room for additional gluonic radiation, and that are color singlets at lowest order (true 
QCD processes can be treated in a similar fashion 0|§). Examples are the e + e~ total 
cross section to hadrons near unit thrust, deep-inelastic scattering near unit Bjorken 
x, Drell-Yan production near threshold, etc. Let us call the edge of phase space 
in such cases the elastic limit. In the integral over virtual and real gluonic degrees 
of freedom, which takes place in the cross section, the most important momentum 
regions are, e.g. in Drell-Yan: (i) fast, almost collinear partons in the two incoming 
jets, (ii) far off-shell, short- distance partons that result from the collision of the two 
incoming jets, and produce the off-shell vector boson, and (iii) soft gluons that couple 
the two incoming jets with momenta p\ and p 2 - Near the elastic limit, the cross section 
factorizes M into corresponding hard, a soft, and two jet functions that summarize 
these degrees of freedom 

a(N) = / dwe- Nw a(w)=H( Pl /^p 2 /f,,Ci) S(Q/fiNXi) 
Jo 

xJi(pi ■ Cx/ft, Q/ftN) J 2 (P2 • C 2 /V, Q/ftN) , (3) 

Let us discuss the variables that occur in the above equation. Q is the hard scale 
of the process (the invariant mass of the vector boson produced in Drell-Yan) and 
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w is a dimensionless weight function that is defined to vanish in the elastic limit, is 
insensitive to collinear splittings of partons, and to additional soft radiation and is 
additive near the edge of phase space, i.e. w = Wi+w 2 +w s where Wi are contributions 
to the weight function from momenta in the jets, and w s is the contribution from 
soft momenta. In DY one can choose w = 1 — z = 1 — Q 2 /s. Further, N is the 
moment variable Laplace-conjugate to w, and is large near the elastic limit. Arbitrary 
variables, on which the physical cross section may not depend, are most generally the 
factorization scale \x and the vectors Q, which are necessary to define the jet functions 
Jj. The latter can be thought of as gauge-fixing vectors, as one may compute the 
jet functions in different gauges, as long as the total cross section does not depend 
on these gauges. It is the arbitrariness in these vectors that allows one to extend 
the earlier arguments for the single logarithm case to the Sudakov double logarithmic 
one. 

One proceeds by acting on this factorized form (|3|) with two differential operators 
with respect to the arbitrary variables just mentioned: 

A , Pi't-A-p (* = 1,2). (4) 
a/x dpi ■ Q 

Using separation of variables, and solving the (double) differential evolution equations 
obtained with the above operators, one arrives at 0: 

In a(N) = C (0) + (a s /n) In 2 N + In iV + C {1) 1 (5) 

where B^\C^ have to be fitted from a next-to- leading order calculation. As 
expected, the resummed cross section is an exponential, with at most double loga- 
rithms, while higher powers of In N in the exponent arise only from the expansion 
of the running coupling. In the next section we consider a specific and important 
example, viz. the resummation of Sudakov logarithms in Higgs production at the 
LHC. 



2. The Resummed Higgs Production Cross Section 

The presently allowed SM Higgs mass ranges from about 70 GeV (from direct 
searches at LEP 0) to about 700 GeV from unitarity/triviality constraints 0. In the 
MSSM the lower limits on the two scalar Higgs bosons h and H, and the pseudoscalar 
A are about 60 GeV ||. The theoretical upper limit on the h mass is about 130 GeV 
|7]]. At the LHC the dominant production process is gluon-gluon to Higgs via a 
top (and bottom, about 10%) quark loop. Two approximations are in order. First, 
we shall neglect initial states involving quarks (they contribute only 10% at NLO). 
Second, we would like to consider the gluon-gluon-Higgs coupling as effectively 
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pointlike. This can be done by taking the infinite top mass limit, supplemented by 
low-energy theorems ||. Comparing the NLO infinite top mass limit result |§ with 
the full analytic massive NLO result [|IIJ] one finds a difference less than 10 % for the 
full range for the SM, as well as for the MSSM provided tg/3 is not too large. We may 
then write for the ci-dimensional partonic cross sections (d — 4 — 2e) 



a 



yy 



(6) 



with the coefficients 



On 



h,H 



9t 



h G F a 2 sB N c C F T 2 (l + e) /4tt 



2c 



115272- 



7T 



(7) 
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A G F al B N c C F T 2 (l + e) 



9t 
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Ait 



2c 



512v^tt 

where a Sy B is the bare strong coupling constant (with dimension 2e) and gf{4> = 
h, H, A) denote the modified top Yukawa couplings normalized to the SM coupling, 
which are given in ||10|| . The effective coupling constants are given to NNLO in 



TT|.[T2|. The correction factor may be expanded perturbatively 

oo 

p^,M|//i 2 ,a(/i 2 ),e) = E cr(f)pf{z,M$ff,e) 



(9) 



n=0 



where we define a(p 2 ) = a s (p 2 )/TT. The lowest and next order components of p are 

= 6(1 -z) 



(10) 



P^Ml/p^e) 
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Pa 



A + 2 [jjs) C A 5{l-z) 



// 



We have implicitly redefined the scale \x by p 2 — > p 2 exp[— (ln(47r) 



{1-zf^} (11) 
(12) 

7e)], and scaled 



an overall 1/z into the parton distribution functions JTT 



Let us now construct a resummed expression for p<j,(z, M?//i 2 , a, e), or rather for its 



Mellin tranform p<f>(N, 



Jq z n 1 p<j,(z, . . .). We could follow the methods outlined 



in the first section, but can in fact simplify further. Assuming the Higgs cross section 
factorizes near the elastic edge of phase space, one would arrive at a resummed cross 
section of the form (||). By acting on this equation with (M^d/dM 2 ) one trivially 
obtains an evolution equation in the Higgs mass of the form: 
d 



Mi 



dM\ 



r^(7V, M£/ (i , a(p 2 ) } e) = W^(N } M 2 /p 2 , a(p 2 ) } e)~p <p {N, M 2 /p 2 , «(// 



(13) 
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In order to solve eq. ( |13D we must impose a boundary condition, and find the evo- 
lution kernel. We may in fact use PJTlH the boundary condition p</>(N, Ml/p 2 = 



0, a(/i 2 ), e) = 1 for the moments, or in z-space p<f,(z, M 2 /p 2 = 0, a(/i 2 ), e) = 6(1 — z). 
The solution to eq. (|l|) is then 



(N,M 2 /p 2 ,a(p 2 ),e) = exp 



which may be expanded, using renormalization group invariance, as 



(14) 



fr(N, ^a(p 2 ),e) = exp[ ^ ^ dX 

p 2 



/q dzz N - l j Q > ^{a(\,a(p 2 ),e)w£\z,l/v,e) + ...} 

(15) 

with v(z) an arbitrary function. The one loop coefficient of the evolution kernel 
can then be derived from a low order calculation of the correction factors via 

W^(N,l,e) = (M 2 /p 2 YM 2 ^p { ;\N,Ml/p 2 ,e) , (16) 



The general structure of the result is, in z-space 



W$\z, 1, e) = 6(1 - z)ff\e) + [ { l_ Z z yi e ) +h {1) (z,e), (17) 

where the coefficient functions /i , g^, are regular functions of their arguments 
at z — 1. We drop the term as it is of order 1/iV 4 in moment space. After 
rescaling to incorporate the factor (1 — z)~ 2e and combining the plus distribution with 
the Mellin transform in eq. fllB"!) , we see that the relevant function to approximate is 
(z N ~ l — l)gW(z,e). We approximate the residue function g^(z,e) in three schemes, 
which are defined by 



scheme a 


1 


[z N ~ l 


-l)9 (1) (z,e) 


- - 1) 
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scheme f3 


1 

c~ A 




-l)gW(z,e) 


- - 1) 
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(l-z)(2z 2 


-Az- 2z 3 


scheme 7 


1 

C~ A 


[z N - 1 


-l)gW{z,e) 


- - 1) 


2 - 


(l-z)(2z 2 


-Az- 2z 3 



-Az N -\l-z). (18) 

The minimal scheme a involves replacing g^(z,e) simply by g^(l,e), scheme (3 
includes all terms of (9(1) in the exponent, whereas scheme 7 includes in addition some 
0(\n l N/N) terms. Using the one- loop evolution kernel we can now construct 
the resummed expressions for the Higgs production correction factor in the three 
schemes. Although these expressions are still divergent for e — > 0, the divergences 
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may be cancelled by mass factorization and renormalization for which we chose the 
MS scheme. 

The final results for the resummed cross sections in moment space are given in 



nj, both for Higgs production and the Drell-Yan process. These two processes are 
very similar from the soft gluon point of view, the main difference being that Higgs 
production is driven by gluon fusion, and Drell-Yan by quark-antiquark annihilation. 
Rather than evaluate the resummed answers numerically - which involves the difficult 
problem of treating the infrared renormalon - we expanded the resummed answers 
|TT| for Higgs production and Drell-Yan |T3| to NNLO in QCD perturbation theory, 
in the above three schemes. The answers are expressed in terms of plus-distributions 
T>i{z) and logarithms £i(z) (which are integrable but large near the edge of phase 
space) 

= [ ln y_~ Z) 1 , Si(z) = ln*(l - z), (19) 

and constants. Scheme 7 incorporates the £ logarithms, which have not been included 
in resummed cross sections before. At any order, the leading £ logarithms and those 
subleading ones that are related to the QCD running coupling are universal, as they 
arise directly from the splitting function. They occur in exact NNLO calculations for 
DY and DIS Jl4|], and we checked analytically that the scheme 7 resummed answer 



for these processes, expanded to NNLO, reproduces them. In Fig. [I] we present 
the correction factors for SM Higgs production at the LHC, which coincide with 
the correction factors of MSSM scalar Higgs boson production for small tg/3. For 
MSSM pseudoscalar Higgs production we show similar results in [II|. In Fig. [TJa 



the "partonic" K-factors, obtained from folding the correction factors with NLO 
parton densities and using a NLO strong coupling for all orders of the cross sections, 
are presented. For comparison we show in Fig. |I|b the corresponding NLO "hadronic" 
K-factors normalized to the LO cross sections evolved with LO parton densities and 
a s . Whereas the former indicate the rate of convergence of the individual-order 
contributions within a fixed order calculation, the latter exhibit the convergence of 
the perturbative approach to the physical (hadronic) quantities. We observe from 
Fig. [I] that at NLO scheme 7, remarkably, reproduces the exact NLO calculation 
almost exactly for the full range of the SM Higgs mass Mjj ?t 70 GeV (similar results 
are obtained for the neutral Higgses in the MSSM) schemes a and f3 agree with the 
exact result only for M ^> 1 TeV (the agreement of scheme a in the intermediate 
Higgs mass range is accidental). Moreover, note that the NNLO corrections to the 



partonic cross sections in scheme 7 are still very significant. We checked |Tl] in Drell- 



Yan that similarly good agreement is obtained in scheme 7 in NLO and NNLO for 



the same kinematics [|14|| . Keep in mind that full NNLO predictions for hadronic 
cross sections require NNLO parton densities, which are not yet available. Should 
the size of the NNLO corrections to the physical cross section warrant concern about 
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Figure 1: a) Exact (solid line) and approximate two- (a\,(3i,ji) and three-loop 
(0.2,(^2,12) partonic K-factors, convoluted with the NLO gluon-gluon luminosity 
dC 9 ^ LO /dr, where t = M%/s, in the heavy top-mass limit and in three different 
schemes, versus the scalar Higgs mass Mh- We used NLO CTEQ4M parton den- 
sities £3} and a s (A|| = 202 Me V). b) Hadromc NLO K-factor using LO CTEQ4L 
parton densities |7^/ and a s (Ajfo = 181 MeV) for the LO cross section and including 
the NLO contributions from kh. 

the convergence of the perturbative approach, our resummed answer provides a tool 
to control such large corrections. 
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